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Abstract. The +49A/G polymorphism of the cytotoxic 
T-lymphocyte-associated antigen-4 gene (CTLA-4) has been 
associated with Graves' disease (GD). However, results have 
been inconsistent. The aim of this study was to quantitatively 
summarize the evidence for CTLA-4 +49A/G polymorphism 
and GD. Electronic search of PubMed was conducted to select 
studies. Case-control studies containing available genotype 
frequencies of CTLA-4 +49 were chosen, and odds ratio (OR) 
with 95% confidence interval (CI) was used to assess the strength 
of this association. Forty-two case-control studies including 
8,288 cases and 9,372 controls were identified. Three studies 
were eliminated from the total 42 studies due to a p-value <0.05 
(p-value for Hardy-Weinberg equilibrium in control group) in 
these studies which induced significant publication bias. The 
overall results suggested that the variant genotypes were highly 
associated (p<0.01) with GD risk in all genetic models (addi- 
tive model: OR, 1.443; 95% CI, 1.319-1.578; p<0.001; recessive 
model: OR, 1.589; 95% CI, 1.396-1.808; p<0.001; dominant 
model: OR, 1.621; 95% CI, 1.430-1.837; p<0.001). Similarly, 
in the subgroup analyses for ethnicity (Caucasian, Asian), 
the results were positive. This meta-analysis suggests that the 
CTLA-4 +49A/G polymorphism is highly associated (p<0.01) 
with increased risk of GD, especially in Caucasians and Asians. 
To validate this association, further studies with larger partici- 
pants worldwide are needed to examine associations between 
this polymorphism and GD. 

Introduction 

Graves' disease (GD) is one of the autoimmune thyroid diseases 
(AITDs) which affect 5% of the general population (1). GD is an 
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autoimmune antibody-mediated, thyroid-specific autoimmune 
disease which causes thyroid gland tumefaction. GD patients 
make antibodies to the thyroid-stimulating hormone receptor 
leading to hyperthyroidism. People of Western countries 
(-1.2%) and 0.25-1.09% of people of China are afflicted with 
GD (2,3). Although environmental factors, such as infection (4) 
and stress, are very important in the process of Graves' disease 
in susceptible individuals, one study in twins revealed that 
~80% of the predisposition to GD is due to genetic factors (5). 
Several genetic loci have been implicated in the susceptibility 
to this disease. One of the associated genes is the cytotoxic 
T-lymphocyte-associated antigen-4 (CTLA-4) gene which 
consists of 4 exons and 3 introns. In 1997, Yanagawa et al (6), 
Marron et al (7) and Donner et al (8) initially reported that 
there was an association between CTLA4 and Graves' disease. 
The CTLA-4 gene is located on the long arm of chromosome 
2q33 and belongs to the immunoglobulin superfamily. Since 
the CTLA-4 protein transmits an inhibitory signal to T-cells, it 
has a strong susceptibility in autoimmunity. One of the CTLA-4 
gene polymorphisms is located on exon 1 +49, which causes a 
threonine to alanine substitution in codon 17 (codon 17 TA). 
To date, the CTLA-4 +49A/G polymorphism has been studied 
in different and numerous groups in humans, and a potential 
association with GD has been found in many results (6-36). 
However, some results suggest that there is no association 
between CTLA-4 +49A/G polymorphism and GD (37-46). 
Thus, the results are still inconsistent. Another problem is that 
these published studies only refer to a rather modest sample size 
that limits their significance. Utilizing the advantage of meta- 
analysis, a powerful method for quantitatively summarizing 
different study results, we combined the data for analysis and 
increased the sample size to a reasonable level. In this study, we 
conducted a meta-analysis to quantitatively assess the effect of 
the CTLA-4 +49A/G polymorphism on the risk of GD. 

Materials and methods 

Publication search. PubMed was searched using the terms 
'CTLA 4', 'Graves' and 'polymorphism' or 'CTLA4', 'Graves' 
and 'polymorphism' or 'cytotoxic T lymphocyte', 'Graves' 
and 'polymorphism' (the last search update was on March 11, 
2012). Case-control studies containing available genotype 
frequencies of 49A/G were chosen. Additional studies were 
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Study 

ID 

Asian 

Kimkonq, L (2011) 

Zhao, S.X. (2010) 

Kimura. H. (2009) 

Wang, P. W. (2007) 

Zhang. Q. (2006) 

Han, §. Z. (2006) 

Cho, H. J. (2006) 

Weng. Y C. (2005) 

Iwama, S. (2005) 

Wang. P. W. (2004) 

Mochizuki. M. (2003) 

Bednarczuk. T. (2003) 

Km w, Y. (2002) 

Parte. Y. i. (20W) 

Awata.T. (1998) 

Yanagawa. T. (1997) 

Marron, M. P. (1997) 

Subtotal (l-squared = 55.5%, p = 0.003) 

Causasian 

Khalilzadeh. 0. (2009) 
Esteghamati. A. (2009) 
Tanrikulu, S. (2006) 
Sahin. M. (2005) 
Petrone. A. (2005) 
Nakkash-Chmaisse, H. (2004) 
Frydecka. I. (2004) 
Vaidya. B. (2003) 
Kalantari. T. (2003) 
Bednarczuk. T (2003) 
Kouki, T. (2002) 
HadjKacem, H. (2001) 
AOahabadia. A. (2001) 
Chistyakov, D. A. (2000) 
Heward. J. M. (1999) 
DjilaS-Saiah. I. (1998) 
Braun, J. (1998) 
Donner, H. (1997) 

Subtotal (l-squared = 74.0%, p = 0.000) 
Others 

Namo Cury, A. (2008) 
Villanueva. R. (2000) 
Kouki, T. (2000) 
Chen, Q. Y. (2000) 

Subtotal (l-squared = 0.0%, p = 0.845) 



Overal (l-squared = 66.1%, p = 0.000) 

NOTE: Weights are from random effects analysis 
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Figure 1. Forest plot of ORs of the G allele when compared to the A allele (additive model) in the Graves' patients. The squares and horizontal lines correspond 
to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight. The diamond represents the pooled OR and 95% CI. OR, odds 
ratio; CI, confidence interval. 



identified by a manual search of the references of the original 
studies. 

Statistic analysis. For the control group of each study, the 
observed genotype frequencies of the CTLA-4 +49A/G 
polymorphism were assessed for Hardy-Weinberg equilib- 
rium using the x 2 test. The strength of association between 
the +49A/G polymorphism of the CTLA-4 gene and GD 
was assessed by calculating crude odds ratios (ORs) with 
95% confidence intervals (CIs). The pooled ORs were 
performed for the additive genetic model (G vs. A), domi- 
nant model (G/G+G/A vs. A/A) and recessive model (G/G 
vs. G/A+A/A), respectively. Heterogeneity assumption was 
checked by a x 2 -based Q-test. A p-value of <0.05 for the 
Q-test indicated a lack of heterogeneity among the studies; 
the summary OR estimate of each study was calculated by 
the random effects model (47,48). The potential for publica- 
tion bias was examined by Begg's test (funnel plot method) 
and Egger's linear regression test (p<0.05 was considered 
representative of statistical significance) (49). All statistical 



analyses were performed with Stata software (version 11.0; 
Stata Corporation, College Station, TX). 

Results 

Eligible studies. We identified 42 case-control studies 
concerning the association between the CTLA-4 +49A/G poly- 
morphism and GD, which included 8,288 GD cases and 9,372 
controls. These data were used in our meta-analysis (Table I). 
The distribution of genotypes in the controls of all the studies 
was in agreement with Hardy-Weinberg equilibrium. 

Meta-analysis. The results of the association between the 
CTLA-4 +49A/G polymorphism and GD and the heteroge- 
neity test are shown in Table II. The overall results suggest 
that the variant genotypes were highly associated (p<0.01) 
with GD risk in all genetic models [additive model: OR, 1.443; 
95% CI, 1.319-1.578; p<0.001 (Fig. 1); recessive model: OR, 
1.589; 95% CI, 1.396-1.808; p<0.001 (Fig. 2); dominant model: 
OR, 1.621; 95% CI, 1.430-1.837; p<0.001 (Fig. 3)]. Similarly, 
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Table I. Distribution of the CTLA-4 +49A/G genotype for patients with Graves' disease and the controls. 



Population 


Ethnicity or 
origin 


Study 






GD 






Control 










Year 


A/A 


A/G 


G/G 


A/A 


A/G 


G/G 


P-value" 


Caucasian 


South Indian 


Veeramuthumari et al (9) 


2011 


11 


37 


32 


29 


25 


26 


0.000819 


Asian 


Thai 


Kimkong et al (37) 


2011 


22 


49 


61 


26 


73 


54 


0.875319 


Asian 


Chinese Han 


Zhao et al (10) 


2010 


104 


730 


1030 


156 


823 


945 


0.211832 


Asian 


Japanese 


Kimura et al (11) 


2009 


62 


143 


210 


142 


358 


295 


0.067982 


Caucasian 


Iranian 


Khalilzadeh etal (12) 


2009 


48 


43 


14 


75 


25 


3 


0.606930 


Caucasian 


Iranian 


Esteghamati et al (13) 


2009 


114 


71 


20 


75 


25 


3 


0.606930 


Others 


Brazilian 


Namo Cury et al (38) 


2008 


43 


58 


15 


39 


32 


7 


0.905523 


Asian 


Chinese 


Chong etal (14) 


2008 


7 


73 


97 


24 


56 


71 


0.028090 


Asian 


Taiwanese 


Wang etal{\5) 


2007 


15 


69 


124 


18 


77 


97 


0.633099 


Asian 


Chinese 


Zhang et al (16) 


2006 


2 


29 


58 


7 


26 


27 


0.846451 


Caucasian 


Turkish 


Tanrikulu et al (17) 


2006 


48 


38 


11 


42 


34 


14 


0.120930 


Asian 


Chinese 


Han etal (18) 


2006 


33 


95 


135 


32 


89 


75 


0.520341 


Asian 


Korean 


Cho et al (39) 


2006 


16 


112 


160 


30 


197 


244 


0.240107 


Asian 


Taiwanese 


Weng et al (40) 


2005 


8 


53 


46 


15 


58 


28 


0.091603 


Caucasian 


Turkish 


Sahin et al (19) 


2005 


29 


33 


15 


43 


48 


7 


0.189953 


Caucasian 


Italian 


Petrone et al (20) 


2005 


59 


68 


23 


139 


138 


24 


0.201228 


Asian 


Japanese 


Iwama et al (41) 


2005 


17 


25 


1 


78 


88 


34 


0.287293 


Asian 


Taiwanese 


Wang etal (21) 


2004 


18 


72 


81 


11 


50 


87 


0.316477 


Caucasian 


Lebanese 


Nakkash-Chmaisse et al (22) 


2004 


8 


23 


3 


24 


14 


0 


0.163933 


Caucasian 


Polish 


Frydecka et al (42) 


2004 


32 


50 


17 


50 


84 


20 


0.096480 


Caucasian 


White 


Vaidya et al (23) 


2003 


88 


139 


74 


146 


158 


45 


0.825642 


Asian 


Japanese 


Mochizuki et al (24) 


2003 


1 


6 


13 


12 


27 


21 


0.539129 


Caucasian 


Iranian 


Kalantari et al (25) 


2003 


21 


49 


20 


30 


53 


30 


0.510214 


Caucasian 


Polish 


Bednarczuk et al (26) 


2003 


75 


123 


66 


77 


85 


32 


0.303455 


Asian 


Japanese 


Bednarczuk et al (26) 


2003 


28 


140 


151 


15 


63 


34 


0.093804 


Asian 


Chinese 


Yung et al (27) 


2002 


3 


54 


66 


23 


59 


76 


0.046372 


Caucasian 


USA 


Kouki et al (28) 


2002 


22 


67 


31 


30 


36 


14 


0.576150 


Asian 


Japanese 


Kinjo et al (32) 


2002 


32 


62 


50 


38 


46 


26 


0.107271 


Caucasian 


Tunisian 


Hadj Kacem et al (30) 


2001 


31 


63 


50 


26 


94 


85 


0.998814 


Caucasian 


UK 

Caucasian 


Allahabadia et al (43) 


2001 


136 


262 


86 


192 


198 


34 


0.081624 


Others 


African, 
American, 
Hispanic, 
Asian 


Villanueva et al (44) 


2000 


42 


67 


28 


53 


52 


16 


0.568526 


Asian 


Korean 


Park etal (31) 


2000 


5 


35 


57 


26 


75 


98 


0.061219 


Others 


Not specified 


Kouki et al (32) 


2000 


8 


29 


8 


15 


23 


5 


0.390573 


Caucasian 


Moscow 


Chistyakov et al (33) 


2000 


6 


22 


50 


25 


38 


30 


0.081864 


Others 


African- 
American 


Chen et al (45) 


2000 


20 


25 


4 


23 


19 


5 


0.718804 


Caucasian 


UK 


Heward et al (34) 


1999 


122 


192 


65 


164 


171 


28 


0.067423 


Caucasian 


White 


Djilali-Saiah et al (46) 


1998 


23 


37 


13 


47 


37 


16 


0.069793 


Caucasian 


German, 
Canadian 


Braun et al (35) 


1998 


22 


56 


25 


52 


48 


21 


0.096985 


Asian 


Japanese 


Awata et al (36) 


1998 


11 


44 


57 


58 


197 


170 


0.938310 


Asian 


Japanese 


Yanagawa et al (6) 


1997 


11 


64 


78 


34 


88 


78 


0.287293 


Asian 


Chinese 


Marron et al (7) 


1997 


1 


11 


16 


6 


39 


49 


0.632129 


Caucasian 


German, 
Canadian 


Donner et al (8) 


1997 


81 


161 


63 


135 


149 


41 


0.990935 



"p-value for Hardy- Weinberg equilibrium in the control group. GD, Graves' disease. 
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Table II. ORs and 95% CI for the CTLA-4 +49A/G polymorphism for different genetic models in patients with Graves' disease. 
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<0.001 
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0.965 


Dominant 
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1.431 


; 1.227- 1.670) 


<0.001 


0.349 


0.187 


0.196 


(G carriers vs. A/ A) 
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1.727 
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OR (95% Cl) 
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Zhao, S X. (2010) 
Kimura. H. (2009) 
Wang, P. W. (2007) 
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Han, S. Z. (2006) 
Cho, H. J. (2006) 
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Khaiilzadeh, 0. (2009) 
Esteghamati. A. (2009) 
Tanrikukj, S. (2006) 
Sahin. M. (2005) 
Petrone, A. (2005) 
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Namo Cury, A. (2008) 
Viilanueva. R. (2000) 
Kouki, T. (2000J 

Chen, Q. V. (2000) 

Subtotal (l-squared = 0.0%, p = 0.773) 

Overall (l-squared = 57.4%, p = 0.000) 
NOTE: Weights are from random effects analysis 
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Figure 2. Forest plot of ORs of the G/G genotype when compared to the A allele carriers (GA+AA) (recessive model) in the Graves' patients. The squares and 
horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight. The diamond represents the pooled 
OR and 95% CI. OR, odds ratio; CI, confidence interval. 



in subgroup analyses for ethnicity (Caucasians, Asians), the Publication bias. Funnel plot and Egger's test were performed 
results were positive. to estimate the publication bias of studies. The results of Egger's 
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Study 
ID 



OR (95% CI) 



% 

Weight 



Asian 

Kimkong, 1.(2011) 
Zhao, S. X (2010) 
Kimura, H. (2009) 
Wang, P. W. (2007) 
Zhang, Q. (2006) 
Han, 5. Z. (2006) 
Cho, H. J, (2006) 
Weng, Y. C. [2005) 
Iwama, S. (2005) - 

Wang, P. W. (2004) ■ 

Mochizuki, M. (2003) 
Bednarczuk. T. (2003) 
Kinjo, Y. (2002) 
Park. Y.J. (2000) 
Awata, T. (1998) 
Yanagawa. T. (1997) 

Matron, M. P. (1997) 

Subtotal (l-squared ■ 9.0%, p ■ 0.349) 

Causasian 

Khalzadeh. 0. (2009) 
Esteghamati, A. (2009) 
Tanrikulu. S. (2006) 
Sahin, M. (2005) 
Petrone, A. (2005) 
Nakkash-Chmatsse. H. (2004) 
Frydecka. I (2004) 
Vaidya, B. (2003) 
Kalantari, T. (2003) 
Bednarczuk, T. (2003) 

Kouki, T. (2002) _ 
HadjKacem, H. (2001) — * 
Alahabadia, A. (2001) 
Chistyakov. D. A (2000) 
Heward.J. M. (1999) 
Djilalt-Saiah, L (1998) 
Braun, J. (1998) 
Donner, H. (1997) 

Subtotal (l-squared = 65.9%, p = 0.000) 
Others 

Namo Cury. A. (2008) 
VBanueva. R. (2000) 
Kouki, T. (2000) 
Chen, Q. V. (2000) 

Subtotal (l-squared ■ 0.0%. p - 0.850) 
Overal (Wsquared = 48.0%. p = 0.001 ) 

NOTE Weights 3re fr;-n tj-scti effsas a-alvsij 




o 

o 



1.02(0.65. 
1.49(1.16. 
1.24(0.89. 
1.33(0.06. 
5.75(1.15. 
1.36(0.80. 

1 16(0.62. 

2 16(0.87. 
0.98 (0.50. 
0.68(0.31. 
4.75 (0.58. 
1.61 (0.82. 
1.85 (1.06. 
2.77 (1.03. 
1.45(0.73. 
2.64 (1.29. 
1.84(0.21. 
1.43(1.23. 



3.18(1.78. 
2.14(128. 
0.89 (0.50. 
1.29(0.70. 
1.32(0.89. 
5.57(1.99. 
1.01 (0.59. 
1.74 (1.26. 
1.19(0.62. 

1.66 (1.12. 

2.67 (1.40. 
0 53 (0.30. 
2.12(1.61. 
4.41 (1.71. 
1.74 (1.29. 
1.93 (1.03. 
2.77 (1.53. 
1.96 (1.40. 
1.73(1.42. 



1.91) 

1.93) 

1.71) 

2.72) 

28109) 

2.30) 

2.16) 

5.34) 

1.92) 

1.50) 

39.10) 

3.13) 

3.22) 

7.44) 

2.87) 

5.41) 

15.97) 

1.67) 



f.631 
3.63'; 
1.59) 
2.38) 
1.97) 
15.62) 

I. 73) 
2.41) 
2.26) 
2.46) 
5.10) 
0.94) 
2.79) 

II. 42) 
2.34) 
3.62) 
5.02) 
2.75) 
2.10) 



1.70(0.95. 3.04) 
1.76(1.06. 2.94) 
2.48 (0.92. 6.66) 
1.39(0.62. 3.12) 
1.74(1.25.2.41) 



1.62(1.43.1.84) 100.00 



2.4? 

4.se 

4.40 
2.06 
0.55 
2.98 
2 44 
1.46 
2.22 
1.81 
0.33 
2.25 
2.80 
1.27 
2.19 
2.05 
0.32 
36.68 



2.67 

1. y. 

2.70 
2.52 
3.83 
1.19 
2.90 
4.39 
2.36 
3.SS 
2.34 

2.7: 

4S- 

■ V 
4.62 
2.42 

2. eo 

4.31 
54.68 



2.66 
3.07 
1.28 
1.74 
8.74 



1~ 
10 



-1 



Figure 3. Forest plot of ORs of GD G allele carriers (G/G+G/A) when compared to the A/A genotype (dominant model) in the Graves' patients. The squares 
and horizontal lines correspond to the study-specific OR and 95% CI. The area of the squares reflects the study-specific weight. The diamond represents the 
pooled OR and 95% CI. OR, odds ratio; CI, confidence interval. 



test provided statistical evidence for funnel plot symmetry (for 
G/G+G/A vs. A/A, p=0.166) (Table II). 

Discussion 

This meta-analysis examined the association of the 
CTLA-4 +49A/G polymorphism with GD and included 
8,288 GD cases and 9,372 controls. Three studies were 
eliminated from the total 42 studies due to a p-value of <0.05 
(p-value for Hardy-Weinberg equilibrium in control group) in 
these studies which induced significant publication bias. The 
results of Egger's test provided statistical evidence for funnel 
plot symmetry (for G/G+G/A vs. A/A, p=0.166). The overall 
results suggest that the variant genotypes were highly associated 
(p<0.01) with GD risk in all genetic models (additive model: 
OR, 1.443; 95% CI, 1.319-1.578; p<0.001; recessive model: OR, 
1.589; 95% CI, 1.396-1.808; p<0.001; dominant model: OR, 
1.621; 95% CI, 1.430-1.837; p<0.001). Similarly, in subgroup 
analyses for ethnicity (Caucasians, Asians), the results were 
positive. 

GD is a disease with significant clinical consequences. 
The mechanism of GD is still relatively unknown. Although 



environmental factors, such as infection (4) and stress, are 
important in the process of Graves' disease in susceptible 
individuals, one study in twins suggests that ~80% of the 
predisposition to GD is due to genetic factors (5). Single 
nucleotide polymorphisms (SNPs) can be used as a tool for 
investigating genetic variations and disease susceptibility. 
GD is an autoimmune antibody-mediated, thyroid-specific 
autoimmune disease. The CTLA-4 protein can transmit an 
inhibitory signal to T-cells and has a strong susceptibility in 
autoimmunity. CTLA-4 protein has recently been described 
as a gatekeeper of conjugation timing and reduced conjugation 
may protect against prolonged contact periods of cytotoxic T 
lymphocytes with autoantigen-defined targets (50). It has been 
in the centre of attention for its key role in autoimmunity. The 
+49A/G polymorphism is one of the CTLA-4 three forms 
of polymorphisms. To date, a multitude of different studies 
were carried out concerning the association between the 
CTLA-4 +49A/G polymorphism and GD, but the results are 
inconsistent. In many studies (6-36) the results are positive, 
however in others (37-46) the results are negative. 

This meta-analysis revealed a highly significant (p<0.01) 
association between the CTLA-4 +49A/G polymorphism and 



Sletal: CTLA-4 AND GRAVES' DISEASE 



543 



GD risk, in both Asian and Caucasian subgroups. In conclu- 
sion, this meta-analysis suggests that the CTLA-4 +49A/G 
polymorphism is potentially associated with the risk of GD 
among Caucasians and Asians. Future, well-designed, large 
scale studies are necessary to validate this association in 
different populations. 
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